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ABSTRACT 

We present our analysis of the red-channel CoRoT data of extrasolar planet CoRoT-2b. A deep secondary eclipse is detected at a level 
of 1.02 + 0.20 X lO""*, which suggests that all of the planet-signal detected previously in white light by Alonso et al. (2009) originates 
from the red channel. CoRoT-2b is the coolest exoplanet that has been detected in the optical so far. In contrast to the other planets, 
its measured brightness temperature of 2170 + 55 K is significantly higher than its maximum hemisphere-averaged effective day-side 
temperature. However, it is not expected that a hot Jupiter radiates as a black body, and its thermal spectrum can deviate significantly 
from a Planck curve. We present models of the planet/star flux ratio as function of wavelength, which are calculated for a T/P profile 
in radiative and hydrostatic equilibrium, using a self-consistent atmosphere code. These are compared with the CoRoT detection, 
and with measurements at 4.5 and 8 fim from the Spitzer Space Telescope from Gillon et al. (2009). We estimate that reflected light 
contributes only at a 10-20% level to the total optical eclipse depth. The models allow for an 'extra absorber' to be inserted at high 
altitude in the planet's atmosphere. This produces a thermal inversion layer, recently invoked to explain the photometric reversals and 
flux enhancements seen in some planets in the infrared. In the 0.5-1.5 yum wavelength range, the model-spectra of planets with an 
extra absorber at high altitude, are relatively suppressed in flux compared to those witliout such absorber. We therefore argue that, 
in contrast to the other exoplanets detected in the optical so far, CoRoT-2b may not exhibit a significant thermal inversion in its 
atmosphere, causing its optical brightness temperature to be boosted above its maximum effective day-side temperature. 
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' 1. Introduction 

Hot Jupiters, a class of extrasolar planets that orbit their host 
' star at only a few percent of the Earth-Sun distance, are prin- 
, ciple targets for studying optical light, which can be a combi- 
' nation of reflected star-light and/or thermal emission. So far the 
optical reflectivity of these strongly irradiated planets have been 
found to be very low, with stringent upper limits to their geo- 
metric albedo of Ag <0.12 for HD75289b, and Ag <0.12 for 
HD209458b, compared to Ag = 0.41 - 0.52 for the solar sys- 
tem gas giants (Cox 2000). This is as expected from theoretical 
predictions, in particular in the absence of clouds in the planet 
atmospheres. Recently, optical emission from three extrasolar 
planets has been measured, by detection of the planetary eclipse 
of OGLE-TR-56b (Sing & Lopez-Morales 2009), and by detec- 
tion of the phase variations of CoRoT-lb (Snellen, de Mooij, 
& Albrecht 2009, see also Alonso et al. 2009), and HATP-7b 
(Borucki et al 2009). The detected optical radiation from these 
planets is most likely from the short-wavelength tail of their ther- 
mal emission, and not reflected star-light, which shifts into the 
optical wavelength regime for these very hot planets. 

In this paper we report on our analysis of the red chan- 
nel CoRoT data on the secondary eclipse of extrasolar planet 
CoRoT-2b, and investigate whether the deep eclipse is due to 
only thermal radiation or also reflected star-light. CoRoT-2b is 
an object with 3.3 Jupiter masses and 1.465 Jupiter radii, or- 
biting a G dwarf star (Teff = 5625 K) with an orbital period of 
1.743 days (Alonso et al. 2009). Notably, the host star shows 
a profound low frequency modulation in brightness at a level 
of 1-2%, due to the presence of spots on its surface, exhibiting 
a rotational period of 4.52 days, indicating it is a young sys- 
tem possibly <500 Myr old (Lanza et al. 2009). In Section 2 the 



CoRoT observations, data reduction and analysis are presented. 
In Section 3 the red channel eclipse measurement and the 4 and 
8 yU Spitzer observations (Gillon et al. 2009) are compared with 
models of the planet/star flux ratio as function of wavelength, 
produced with a self-consistent atmosphere code and solution 
techniques as described earlier in Burrows et al. (2008). The re- 
sults are discussed in Section 4. 

2. Observations and analyses 

2.1. CoRoT observations 

The first long run of the CoRoT satellite (Auvergne et al. 2009), 
from May until October 2007, was pointed towards the galactic 
center direction, monitoring the flux of ~12,000 stars. CoRoT-2b 
was identified as a planetary candidate by the alarm mode within 
a few days (Alonso et al. 2009). This alarm mode is aimed at 
identifying a reduced number of stars for which the sampling 
rate is changed from 512 s to 32 s. In this way, CoRoT-2b was 
observed almost continuously for 142 days, with a sampling rate 
of 5 12 sec for the first 5 days, and 32 sec for the remainder of the 
time, resulting in ~369,700 flux measurements. For our analysis 
we used the updated N2 pipeline datcQ which was released to 
the general public on February 16, 2009. 

A prism in front of the exoplanet-CCDs produces a small 
spectrum for each star, on which aperture photometry is per- 
formed in three bands (red, green, and blue; Auvergne et al. 
2009) This is done on board to comply with the available teleme- 
try volume. The transmission curves of the three bands are differ- 
ent for each targeted star, and depend on the template chosen for 

' http://idc-corotn2-public.ias.u-psud.fr/ 
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Table 1. The star and planet parameters from Alonso et al. 
(2008), which are used for our analysis. Parameter P is the or- 
bital period, T^ is the transit time, T, is the effective temperature 
of the star, M, and R« are the stellar mass and radius, Ry,/R, is 
the ratio of planet to star radius, e is the eccentricity, and I is the 
orbital inclination. 



P[d] 


= 1.7429964+0.0000017 


Tc[d] 


= 2454237.53562±0.00014 


T. [K] 


= 5625+120 


M. [M,„„] 


= 0.97+0.06 


R* [Rsun] 


= 0.902+0.018 


Rp/R. 


= 0.1667+0.0006 


e 


= (fixed) 


I [deg] 


= 87.84+0.10 



the on-board aperture photometry, which in itself is based on the 
effective temperature of the star and its position on the CCD. For 
CoRoT-2 template ID number 25 was used with the red chan- 
nel covered by 38 pixels. The background was measured on 400 
windows of 10x10 pixels in the field, of which 300 and 100 with 
samphng of 512 s and 32 s sampUng respectively (Auvergne et 
al. 2009). 

To determine the short wavelength cutoff of the red channel, 
we first compared the Kurucz model spectrum (Kurucz 1993) 
for the host star with its apparent magnitude in B, V, r, I, J, H, 
and K bands as taken from the EXODAT and 2MASS databases 
(Meunier et al. 2007; Skrutskie et al. 2006). We found that the 
optical to near-infrared colors are best fitted with an interstellar 
extinction of Av=0.7. A faint star at 4 arcseconds from CoRoT-2 
falls completely within its mask template. From its broadband 
colors and apparent brightness we determine it to be an M- 
dwarf with an effective surface temperature of T~3500 K, 3.5 
magnitudes fainter than CoRoT-2 in V-band. The short wave- 
length cut-off of the red channel pass-band was subsequently 
determined from the overall transmission curve for the tele- 
scope/CCD combination multiplied with the appropriately red- 
dened Kurucz model spectrum of the host star plus a small con- 
tribution from the M-dwarf, which was compared to the fraction 
of photons collected in the red channel. This results in a wave- 
length cutoff of 580 nm, and an effective wavelength of 720 nm. 

2.2. CoRoT data analysis 

The median red-channel flux of the star is ~503,000 e- per 32 
s exposure, with an average background level of 18,700 e", re- 
sulting in a photon noise limit of 0.00144. To avoid systematic 
effects originating from the beginning or end of the light curve, 
or from the change in samphng, we only used the 32 sec sam- 
pled data from the first to the final transit, covering a total of 
~132 days and 78 planetary orbits. Only those flux measure- 
ments which are assigned status = were included in the analy- 
sis (90% of the remaining data), indicating that they are a valid 
measurement, and avoiding those influenced by an energetic par- 
ticle hit or gUtch, or those acquired during the crossing of the 
satellite over the South Atlantic Anomaly. We further rejected 
3cr outlier points relative to the median-smoothed light curve 
(1.4% of the data). 

The Earth has a significant influence on the photometric per- 
formance of the satellite and introduces relevant perturbations 
on time scales of the satellite orbital period (103 min) and the 
24 hr day, due to ingress and egress of the spacecraft from the 




Fig. 1. The red-channel light curve covering 78 planet orbital pe- 
riods of CoRoT-2b, as used for our analysis. The data have been 
binned to 0.005 in phase. Intrinsic stellar variability, probably 
due to spots, completely dominates variability on time scales > 
1 day (see Lanza et al. 2009). The point-to-point scatter, relative 
to the median- smoothed light curve is 5.5x10"''. 



Earth shadow, variations in the gravity field and magnetic field 
(such as the South Atlantic Anomaly), and changes in the levels 
of thermal and reflected Ught from the Earth. Most of these ef- 
fects were corrected for prior to the data release for the general 
astronomy community (Auvergne et al. 2009), however, resid- 
ual effects are still present in the data. The variations in the light 
curve of CoRoT-2 are dominated by intrinsic stellar variability 
at a level of 1-2%, on time scales of a few days (see Fig. 1). The 
residual 24 h instrumental perturbations, as found in the data of 
CoRoT-1, are therefore indistinguishable from the stellar vari- 
ability and cannot be corrected for, but we do correct for the 
perturbations on a 103 -min time scale. 

We corrected for this 103 min perturbation, which is associ- 
ated with the satellite orbital period, as illustrated in Figure 2. 
First we determined the average flux of CoRoT-2 per 103 min. 
period, and constructed a model of the long-term variability by 
fitting a 7th order polynomial through these points (excluding 
those affected by transits), for blocks of two planet orbital peri- 
ods (corresponding to ~5.4 days). This model was used to nor- 
malize the light curve. We subsequently phase-folded the nor- 
malized Ught-curve of these blocks over 103 minutes and box- 
car smoothed over 5 minutes. These curves were used to correct 
each associated block of data. In this way, the relative standard 
deviation was reduced from 0.0067 to 0.0033, still more than 
a factor of two above the photon noise limit. Afterwards, the 
long-term variability in the data was restored by multiplying the 
light curve with the initial model. Hence, in total we determined 
and corrected for the 103-min perturbation for 39 blocks of data, 
each covering 2 planetary orbital periods. In Figure 3 we show 
the total influence of these corrections on the phase-folded data. 
Assuming an uncertainty of 10% on the amplitude of the correc- 
tions, they contribute to the errors in the final light curve at an 
insignificant level of 2x10"^ per data point. 

2.3. Binning, fitting, and phase-folding 

To reduce the noise per data point and expose possible system- 
atic effects, we averaged and binned the data over regular phase 
steps of = 0.005, corresponding to 12.5 minutes. The rela- 
tive standard deviation of this binned data (~ 15,000 points) rel- 
ative to a median- smoothed light curve is 0.00055 (compared 
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Fig. 2. [top panel] Data covering two planetary orbits, phase- 
folded over the orbital period of the CoRoT Satellite (103 min), 
showing the instrumental jitter at a level of ~0.002. The solid 
line indicates the correction for this effect for this section of the 
data, [lower panel] Small part of the light curve showing the 
same 103-min jitter and its correction (solid line). The dashed 
line is our fit to the long-timescale stellar variability which was 
first removed before the instrumental correction was determined, 
after which it was added again. 



to 0.00031 as expected from the photon noise limit). It exposes 
several positions in the light curve which show sudden jumps 
in flux (Fig. 4). When these occur at a relevant section of the 
planetary phase (0.25<^*<0.75), the data taken during those par- 
ticular planetary orbits were not used for further analysis. We 
removed data from 6 of the 78 orbits in this way. To remove 
the influence of the long-term intrinsic variability of the star, we 
least-squares fitted the data for each planet-orbit with a third or- 
der polynomial between phase O.25<0<O.45 and O.55<0<O.75, 
making sure that the fit is not influence by the possible presence 
of a secondary ecUpse at 0=0.5. The best- fit coefficients were 
subsequently used to correct for the stellar variability, by con- 
structing a polynomial curve for the same range in phase includ- 
ing the possible echpse, 0.25<(;*<0.75. To check for systematic 
effects, we varied the range in phase and order of the polynomial 
to fit the long-term variability (see below). Finally, the data from 
all 72 planet orbits were averaged in phase, and their associated 
errors were estimated from the variance between the different or- 



Fig. 3. The overall correction (solid line) for the 103 min pertur- 
bation to the final binned light-curve (grey squares with 1 error 
bars). Although the corrections have an amplitude at a level of 
0.002 (more than an order of magnitude larger than the measured 
ecUpse depth), they largely cancel out over the ~1900 satellite 
orbits. Assuming an uncertainty of 10% in the amplitude of the 
103-min perturbation signal, it contributes to the errors in the fi- 
nal Ught curve at an insignificant level of 2x10"^ per data point. 
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Fig. 4. Six segments of the light curve, each one planet orbital 
period in length, showing sudden jumps in flux. The data from 
these planet orbits were omitted for further analysis. 



bits, typically -6x10 ^ (compared to 3.7x10 ^ as expected from 
the photon noise limit). 
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Fig. 5. The secondary eclipse of CoRoT-2b. The phase folded 
light-curve is shown from 142 days of CoRoT monitoring in its 
red channel. The data are binned over 0.005 in phase. The inte- 
grated flux from the day- side of the planet relative to the stellar 
flux, as determined from the depth of the echpse, is found to be 
1.02+0.20x10""*. This corresponds to a brightness temperature 
of Tb=2170±55 K, which is significantly higher than the maxi- 
mum possible eff^ective day-side temperature for CoRoT-2b. 



2.4. Chi-squared analysis of the eclipse 

The final phase-folded light curve (Fig. 5) shows a distinct dip 
at 0=0.5, exhibiting the correct width to be the planet's sec- 
ondary eclipse. It was fitted with a three parameter model us- 
ing a chi-squared analysis, with K, the eclipse depth, (po, the 
eclipse timing, and Zq, the normalization outside the eclipse, as 
free parameters. The shape of the secondary eclipse as function 
of phase is calculated from the planet and star parameters (Table 
1) using the algorithm of Mandel and Agol (2002). The analy- 
sis was performed on a grid in 3-dimensional parameter space, 
for which the constraints on a single parameter were determined 
by marginalizing over the others. The errors in the data points 
were increased by 14% to obtain a reduced chi-squared of unity. 
In this way, we determined the two relevant parameters to be 
K= 1.09+0. 19 X lO""* and 0o = 0.495 ± 0.003. Since the devi- 
ation of the eclipse timing from (p=Q.5 is not significant, and 
00=0.5 is the likely solution for such a close-in planet, leaving 
the ecUpse timing to vary freely the chi-squared analysis may 
bias K towards higher values. If the eclipse timing is kept at 
00=0.5, the eclipse depth is measured to be K =1.02+0.20x10"^, 
which we use as our final result. Note that Gillon et al. (2009) 
measure a deviation from a circular orbit using Spitzer, but at a 
smaller level, with ecoia>=-0.0029±0.0006, which corresponds 
to 00 = 0.4982 ± 0.0004, ~lcr away from our measurement. 

2.5. Markov-chain Monte Carlo Simulation 

As an extra check for a degeneracy between the parameters, we 
also performed a Markov-Chain Monte Carlo (MCMC) simu- 



lation (Tegmark et al. 2004), using the Metropolis-Hastings al- 
gorithm (Metropolis et al. 1953; Hastings 1970), which should 
result in similar uncertainties in the model parameters as 
above, since the full 3-dimensional parameter space was already 
mapped. Five independent chains were created, each from a ran- 
dom initial position, consisting of 50,000 points, for which the 
first 10% of each chain were discarded to minimize the efi'ect 
of the initial condition. The Gelman & Rubin statistic (Gehnan 
& Rubin 1992) was calculated for each parameter to check the 
convergence and the consistency between the chains, which is 
a comparison between the intra-chain and inter-chain variance. 
The results were weU within 1% of unity, a sign of good mix- 
ing and convergence. The estimated probability distribution of 
K is shown in the middle panel of Figure 6, with the median of 
the distribution indicated by a sohd fine. The dashed fines en- 
close 68% of the results, with equal probabiUty on either side 
of the median. This corresponds to K = 1.04+0.22 x 10"^. The 
probability distribution of 0o is shown in the top panel of Figure 
6. It corresponds to 0o=O. 495+0. 003. The lower panel of Fig. 6 
displays 0o versus K for the MCMC simulation, showing a ran- 
domly selected subsample of 2x10^ points for plotting purposes. 

2.6. Systematic effects 

We further checked for the influence of systematic effects in 
several ways. First of all related to the correction for the stel- 
lar variability, we determined at what level the choice of the 
fitted phase-range, and the order of the polynomial, influence 
the results. We varied the range in phase between 0.2< < 0.8 
and 0.4< <0.6 in steps of 0.025, and the order of the poly- 
nomial between 2 and 4. The results are shown in Figure 7. It 
shows that neither varying the phase-range nor the order of the 
polynomial used in the fit significantly changes the results. Only 
when a small phase-range and 4th order polynomial is used, does 
the eclipse depth deviate towards higher values. We believe that 
this is due to correlated (instrumental) noise still present at short 
time-scales, which can be seen as low-amplitude 'wiggles' in the 
baseline in Figure 5. If the phase-range, over which the polyno- 
mials are fitted, is too small, the baseline is not fitted correctly 
and the resulting eclipse becomes deeper. However, for a phase- 
range from 0.4 to 0.8, the influence of the instrumental 'wiggles' 
on the baseline averages out. We therefore beUeve that a phase- 
range of 0.5 is well chosen. 

The presents of correlated noise could make us underesti- 
mate our errors. We therefore also tried to characterize the uncer- 
tainty using the "residual permutation" bootstrap method (Winn 
et al. 2009). The best fitting model was subtracted from the data, 
after which the residuals were shifted between 1 and 99 points 
and added back to the model fight curve. These light curves were 
refitted for each individual shift. The resulting distributions of 
parameter values (Figure 8) represent separate estimates of the 
uncertainties in the data. With this analysis we find an uncer- 
tainty of 1.9x10"^ for the eclipse depth and 2.4x10"^ for the 
eclipse timing, hence very similar to what was derived above. 



2.7. Analysis of the Green, Blue, and White light data 

We also performed a similar analysis for the CoRoT blue 
channel, green channel, and combined white light curve 
(blue+green+red charmels), as for the red channel data presented 
in this paper. The resulting phase folded light curves around the 
secondary eclipse are presented in Figure 9. The noise in the 
blue and green channels are larger by a factor ~2.5 compared 
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Fig. 7. The fitted eclipse depth (top panel) and eclipse timing 
(lower panel) for different choices of how to remove the long- 
term stellar variability, varying the order of the polynomial (2nd 
= crosses, 3rd = squares, 4th = triangles), and the total phase- 
range over which the variabihty is fit. The arrows indicate the 
parameters used in our analysis, a phase-range of A<p=0.5 (0.25 < 
<0.75), and a third order polynomial. 



n nl I I I 

0.48 0.49 0.50 0.51 

*0 

Fig. 6. The probability distributions of cpo (top panel) and the 
eclipse depth, K (middle panel) from the MCMC simulations. 
The solid line shows the median value. The dashed lines enclose 
68% of the results with equal probabiUty on either side of the 
median. The lower panel shows 0o versus K from the MCMC 
simulations, showing a randomly selected subsample of 2x10^ 
points for plotting purposes. 



eclipse depth of 6.0±2.0xlO"^ is translated to a brightness tem- 
perature of 1910±100 K. However, for the same eclipse depth, 
with the planet and star parameters as given by Alonso et al. 
(2008), we calculate that it corresponds to a brightness tempera- 
ture of Tb=2130^]|Jq K in the white band. Anyway, since 70% of 
the total white light detected by CoRoT originates from the red 
channel, and the white channel detection is at about a 60% level 
of that in the red channel, it suggests that all planet Ught seen in 
the secondary ecUpse of the white light curve originates (within 
uncertainties) from the red CoRoT channel. 



to that for the red channel data, and show no sign of a sec- 
ondary eclipse. The white channel Ught curve has already been 
presented recently by Alonso et al. (2009) who detect the sec- 
ondary eclipse at a depth of 6.0+2.0x10"^. We confirm this re- 
sult at 5.8+2.0 X 10"^, if we use a small phase-range to fit the 
baseUne. However, when we use a larger phase-range, eg. for 
0.25< (p <0.75, we measure an eclipse depth of 4.5±1.5xl0"^, 
a similar variation as seen in the red light curve, probably due 
to residual correlated noise. Note that in Alonso et al (2009) the 



3. Comparison with atmosphere models 

3.1. A simple black body + reflection model 

The measured eclipse depth in CoRoT's red channel corresponds 
to a planet brightness temperature of 2170±55 K, assuming a 
Kurucz model spectrum for the star (with parameters as given in 
table 1). The error includes the incertainty in the effective tem- 
perature of the star. Although it is generally not expected that a 
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and secondary eclipse depth (lower panel) using the residual per- 
mutation bootstrap method. The smooth curves are Gaussians 
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uncertainties in the eclipse timing and depth as from the chi- 
squared and bootstrap methods. 



hot Jupiter radiates as a black body in this wavelength regime, 
we can in first instance compare this with the expected range 
of effective day-side temperatures for CoRoT-2b. If we assume 
a homogeneous effective temperature over the whole day-side 
hemisphere, and if we assume the planet to be completely non- 
reflective at all wavelength (Bond albedo Ab=0), then we ex- 
pect its hemisphere averaged effective day-side temperature to 
be between 1540 K (for a redistribution factorP,, - 0.5, meaning 
that the planet atmosphere very effectively transports stellar heat 
from the day-side to the night-side), and T=1830 K (for a redis- 
tribution factor P„=0). Assuming the most extreme case, a non- 
reflective planet with a dynamics free atmosphere where the ra- 
diation time scale dominates, the temperatures at the sub-stellar 
point are significantly higher, boosting the hemisphere-averaged 
effective day-side temperature to T,nax=1965 K. Even this max- 
imum effective temperature is significantly lower than the mea- 
sured brightness temperature (at 3cr), which if such planet would 
radiate as a black body, would produce a secondary eclipse with 
a depth of 4x10"^ in CoRoT's red channel. 

In addition to the optical detection presented here, CoRoT-2b 
has recently also been observed in the infrared using the Spitzer 



o 

V 

< 



E 
o 

o 




1400 



1600 



1800 
Temperot 



2200 



2400 



Fig. 10. A simple parametrization of the optical flux from 
CoRoT-2b. The planet's uniform day-side temperature is shown 
versus the albedo in the red channel, for a zeroth-order model 
consisting of a thermal black body spectrum plus a constant geo- 
metric albedo over the red channel pass-band. The three vertical 
dashed lines show the day-side effective temperatures, for a non- 
reflective planet with no redistribution of stellar energy towards 
the night-side (Tpn=o = 1830 K), for a non-reflective planet with 
half of the stellar energy transported to the night-side (Tpn=o.5 
- 1540 K), and for a non-reflective dynamics free atmosphere 
(Tmax=1965 K). The dotted lines show the same as the dashed 
lines, but assuming a wavelength-independent Lambert scatter- 
ing, with the Bond albedo, Ag = 3/2 Ag. The grey band covers 
those values allowed by the data within their Icr uncertainties, 
indicating a geometric albedo of 0.06<Ag<0.2 for this simple 
model. However, note that the spectra of neither self-luminous 
gaseous planets nor highly irradiated planets are expected to 
closely resemble Planck curves, thus this in itself is no proof 
for reflected light. 



Space Telescope, showing eclipse depths of 5.10+0.42 x 10"^ 
and 4.1 + 1.1 x 10"^ at 4.5 and 8 micron respectively (Gillon et 
al. 2009). These measurements correspond to brightness tem- 
peratures of T=1805+70 K and T=1325+180 K at 4.5 and 8.0 
micron respectively, significantly lower than in the optical. 

Keeping the important caveats in mind, that the spectra of 
neither self-luminous gaseous planets (and brown dwarfs) nor 
highly irradiated planets are expected to closely resemble Planck 
curves, we show in Figure [TO] the planet's day-side temperature 
versus the albedo in the red channel for a zeroth-order model, 
consisting of a thermal black body spectrum plus a constant ge- 
ometric albedo over the red channel pass-band. The grey band 
indicates those values allowed by the data within their Icr uncer- 
tainties, implying a geometric albedo in the range 0.06<Ag<0.2. 
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Fig. 9. The phase-folded light curves for the blue channel, green channel, red channel, and combined white data. Although the 
secondary eclipse in the white light curve data is formally detected at a 3cr level, at least in our analysis the fitted eclipse depth 
varies significantly depending on the range in phase over which the stellar variability is fitted and removed. 



3.2. Burrows, Ibgui & Hubeny atmosphere code 

We model the ratio of planet to star flux as function of 
wavelength, using the self-consistent atmosphere code and so- 
lution techniques (Hubeny & Lanz 1995; Burrows, Hubeny 
& Sudarsky 2005; Burrows et al. 2007) as employed for 
HD209458b (Burrows, Ibgui & Hubeny 2008). We produced 
four sets of spectra, which are shown in Figure [TT]as planet/star 
flux ratios (left axis) and planet albedos A^ (right axis) as func- 
tion of wavelength. All models are calculated for an atmosphere 
with a solar metallicity. The green and magenta lines indicate 
the spectra for stellar heat redistribution factors of P„ =0. 1 and 
P„=0.3 respectively. The red and brown lines are for the same 
redistribution factors, but now with an extra absorber distributed 
everywhere in the atmosphere at pressures lower than 0.014 
bars, with an opacity of a-=0.01 cm^g ' in the wavelength range 
0.43 ium<A<l.Q fj.m. Such an 'extra absorber' in the optical has 
been invoked to explain the photometric reversals and flux en- 
hancements possibly caused by a thermal inversion as seen in 
several of the highest irradiated planets, including HD209458b, 
which may also be present in CoRoT-2b. 

The CoRoT and Spitzer (Gillon et al. 2009) measurements 
are superimposed on the model spectra. The horizontal bars 
show the expected model-fluxes for these observations, properly 
integrated over the pass-bands. The lengths of the bars indicate 
the effective width of this band. Note that for all models, Ag 
exceeds unity long-ward of ~ 1 fim. Since in this context Ag is 
strictly defined as the ratio of planet to stellar flux at the same 



wavelength, this does not violate energy conservations, because 
in the infrared the planet flux is dominated by absorption and 
re-emission of the star's optical flux. To guide the eye, but keep- 
ing in mind their limited physical meaning, the flux ratios for 
a hypothetical planet with a black body spectrum for an effec- 
tive temperature of 1780 K (P„=0.1) and 1670 K (P„=0.3), are 
indicated by the black dashed and dot-dashed lines respectively. 

In the models without an extra optical absorber, broad ab- 
sorption bands are present in the infrared, mainly from water, 
methane, and CO, and relatively narrow absorption features can 
be seen in the optical from sodium and potassium. The flux 
ratios strongly decrease from infrared to optical wavelengths 
until Rayleigh scattering off H2 and He becomes dominant at 
<600 nm. The P„=0. 1 model spectrum provides a good fit to the 
CoRoT measurement, giving a planet/star flux ratio of 9.5x10"^ 
in CoRoT's red channel, corresponding to geometric albedo of 
0.152. The models that include an extra absorber provide only 
30-40% of the flux seen in the CoRoT red channel. 



4. Discussion and conclusions 

The measured brightness temperature of 2170+50 K in CoRoT's 
red channel is significantly higher than the maximum possi- 
ble hemisphere-averaged effective day-side temperature for this 
planet. However, as we have stressed before, it is not expected 
that such planet would radiate as a black body in this wavelength 
regime, and therefore does not directly imply that the light is 
dominated by reflected star light. The measured flux corresponds 
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Fig. 11. The ratio of planet to star flux (left axis) and geometric albedo (right axis) are shown as function of wavelength, calculated 
using a self-consistent atmosphere code (Burrows et al. 2008). The green and magenta lines indicate model spectra for stellar heat 
redistribution factors of P„=0.1 and P„=0.3, respectively. The red and brown lines are for the same redistribution factors, but now 
with an extra absorber distributed high up in the atmosphere. Superimposed on the model spectra are the CoRoT measurement and 
Spitzer measurements of Gillon et al. (2009) 

. The horizontal bars show the expected model-fluxes for the CoRoT and Spitzer pass-bands. To guide the eye, but keeping in mind their limited 
physical meaning, the flux ratios for a hypothetical planet with a black body spectrum for an efl'ective temperature of 1780 K (P„=0.1) and 1670 
K (P„=0.3), are indicated by the dashed and dot-dashed lines respectively. 



to a geometric albedo of A^=0. 164+0.032, with A^ strictly de- 
fined as above. 

We compared the optical secondary eclipse depth, and 
Spitzer measurements at 4.5 and 8 fim (Gillon et al. 2009) 
with the Burrows, Ibgui & Hubeny (2009) atmosphere models 
for CoRoT-2b, assuming a range of stellar energy redistribu- 
tion factors, and assuming the presence or not of a strong op- 
tical absorber high up in the planet's atmosphere. It is currently 
thought that the presence or absence of such strong absorbers 
in hot Jupiter atmospheres is connected to the level of irradia- 
tion (Burrows et al. 2005; Fortney et al. 2008). For planets that 
receive less stellar flux, this absorber condenses out of the up- 
per atmosphere, e.g. for HD 189733b, while at higher irradiation 
levels this absorber would cause a thermal inversion layer, as 
thought to be the case for HD209458b. It is often quoted that the 
strong optical absorbers TiO and VO are probably the respon- 
sible species (Fortney et al. 2008). However, this has recently 
been questioned. Using a radiative-convective radiative-transfer 
model and a model of particle settling in the presence of turbu- 
lent and molecular diffusion, it is found that it is unlikely that 
VO plays a role in producing an upper atmosphere thermal in- 



version (Spiegel, Silverio, & Burrows 2009). Furthermore, the 
models imply that a heavy species such as TiO can only persist 
in a planet's upper atmosphere with a strong macroscopic mix- 
ing process (Spiegel et al. 2009). 

Purely from its high level of irradiation, in comparison with 
other planets, it would be expected that CoRoT-2b exhibits a 
temperature inversion, implying that the unspecified extra ab- 
sorber is not condensed out. However, the optical eclipse depth 
is only fitted well by a model atmosphere without such extra ab- 
sorber, while the models with this extra absorber underestimate 
the optical flux by more than a factor two. This favours models 
for CoRoT-2b without a significant thermal inversion layer. 

We can estimate the fraction of reflected light in CoRoT's 
red channel for the set of models from the spectral behaviour at 
<600 nm, where the planet flux is dominated by reflected star 
light. It indicates that for those models without an additional 
high-altitude absorber, the reflected light component should con- 
tribute at a level of 1 -2x 10"^, meaning that reflected light forms 
10-20% of the planet flux is CoRoT's red channel. For those 
models with an extra absorber, the reflected light component is 
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significantly lower, at a level of ~5xlO"^ in CoRoT's red chan- 
nel. 

So far, optical secondary eclipses have been detected for 
three other hot Jupiters, OGLE-TR-56b (Sing & Lopez-Morales 
2009), CoRoT-lb (Snellen et al. 2009; Alonso et al. 2009), and 
HATP-7b (Borucki et al. 2009). It is very interesting to com- 
pare the measured brightness temperatures with the maximum 
possible effective day-side temperatures for these planets, as 
shown in Figure [12] The grey band indicates the area for which 
brightness temperatures fall between the maximum hemisphere- 
averaged effective dayside temperature T^ax and Tp_^=o- It shows 
that CoRoT-2b is not only the coolest planet of the four, it is 
also the only planet for which the optical brightness temper- 
ature exceeds Tmax- We mention again that it is not expected 
that the spectra of irradiated hot Jupiters follow a Planck curve, 
and that therefore the direct comparison of brightness tempera- 
ture with the planet's effective temperature has limited physical 
meaning. However, it is an indication that the atmospheric struc- 
ture of CoRoT-2b is different from the other planets, which we 
believe could be due to the absence of a thermal inversion layer 
in CoRoT-2b. Indeed, Figure[TT]shows that the optical spectra of 
planets with an extra absorber at high altitude, are relatively sup- 
pressed in flux compared to those without such absorber. This 
is because in the models with a weak extra absorber in the at- 
mosphere, as a result of such absorption the optical flux can 
be lower than for models without such absorber (Burrows et al. 
2008). As one first raises the amount of extra optical absorber at 
low pressures (high altitudes), the temperature in the upper at- 
mosphere (at the same low pressures) goes up. This is the source 
of the extra mid-infrared flux. At the same time, it absorbs in the 
optical thereby suppressing emergent optical flux. Furthermore 
the temperature at higher pressures goes down, which is the ori- 
gin of some of the near-infrared and also emergent optical flux. 
As one raises the amount of extra absorber further, for large stel- 
lar fluxes at the sub-stellar point, the temperature at low pressure 
continues to go up, as does the mid-infrared flux. Now, the tem- 
perature in the near-infrared photospheric regions also starts to 
rise, which at this point also leads to an increase again of the 
emergent optical flux. Hence, in this regime the optical flux is 
not a monotonic function of the temperature at low pressure, but 
first goes down and then up (Burrows et al. 2008). 
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